Generation of valley polarized current in bilayer graphene 
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We propose a device for the generation of valley polarized electronic current in bilayer graphene. 
By analyzing the response of this material to intense terahertz frequency light in the presence of 
a transverse electric field we demonstrate that dynamical states are induced in the gapped energy 
region, and if the system parameters are properly tuned, these states exist only in one valley. The 
valley polarized states can then be used to filter an arbitrary electron current, so generating a valley 
polarized current. 



One particularly interesting feature of mono- and bi- 
layer graphene [1-3] is the valley degree of freedom. The 
six corners of the Brillouin zone (the K points in the inset 
to Fig. 1) are separated from each other in momentum 
space, and the geometry of the reciprocal lattice requires 
that opposite corners are inequivalent so that there are 
two species of K point, called 'valleys' [4]. The low energy 
spectrum is localized near the six K points, so that in this 
limit, which of the two valleys the electron momentum is 
located in becomes a good quantum number. The valley 
degree of freedom therefore constitutes a two state sys- 
tem (analogous to the electron spin) and is often called 
the 'isospin'. This has prompted the suggestion that the 
isospin could be manipulated and controlled in a useful 
way (so-called 'valleytronics'), for example, to make a 
solid state qubit [5]. Of course, in order to achieve this 
goal, one must be able to accurately prepare and ma- 
nipulate electron states in one valley or another, and to 
date there have been several proposals for devices which 
purport to achieve this [5-10]. 

Recently, attention has also turned to the optical prop- 
erties of monolayer graphene, and its response to lin- 
early and circularly polarized irradiating light fields has 
shown interesting features resulting from the chirality of 
the electrons and the linear low energy spectrum [11, 12]. 

In this Letter, we combine these areas of interest and 
analyze the response of the energy spectrum of gapped 
bilayer graphene [2, 13] to external electromagnetic radi- 
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FIG. 1: Schematic of the valley filtering device. The area 
under the transparent top gate is irradiated, and the parts 
of the flake lying outside of the gated region function as the 
graphitic leads. The inset shows the first Brillouin zone. 



ation in the terahertz frequency range. We then propose 
a device which filters electrons according to which valley 
they are in, creating a valley polarized current. Specifi- 
cally, we find that the different sublattice composition of 
the wave functions of electrons in opposite valleys causes 
them to interact with the irradiating field asymmetri- 
cally. When the radiation and system parameters are 
properly tuned, dynamical states existing entirely in one 
valley are induced. If a current of electrons in this en- 
ergy range is passed through the irradiated region, the 
absence of available states in one valley means that those 
electrons are unable to pass, while electrons in the other 
valley may. The current exiting the irradiated region is 
therefore comprised of electrons in only one valley, a so- 
called 'valley polarized current'. 

This filtering effect is a direct result of the valley asym- 
metric density of states in the irradiated region, and is 
therefore a bulk effect, independent of the geometry of 
the sample and its edges. This gives our device a signifi- 
cant advantage over many prior proposals as it does not 
rely on the precise construction of an edge (as in Refs. 
5-7), or the exact deposition of a gate along one crys- 
tallographic direction (as in Ref. 8), both of which are 
very challenging tasks. Reference 10 also necessitates a 
complex gating arrangement to support one-dimensional 
channels in the graphene. Even if these devices could be 
manufactured, the currents they produce are often only 
partially polarized, and are localized in one-dimensional 
channels, whereas our proposal shows complete valley po- 
larization for significant current flow in a bulk situation, 
making the potential for applications of the current gen- 
erated by this device much more plausible. 

We model irradiated bilayer graphene using the Hamil- 
tonian H = Hq + Hu + H{t), where Hq is the Hamil- 
tonian of ungated, unirradiated graphene and Hjj repre- 
sents the inter-layer potential difference generated by the 
top gate [4]. The time dependent term H(t) is the Hamil- 
tonian of the irradiating field, described by making the 
Peierls substitution in Hq + Hjj with the vector potential 
A = Fy|£!| [cos ilt, sin fit] (where ft is the frequency of 
the radiation) giving 



Hit) 



£,v F eF fQ 1 

mi U o 



e 

3 zf2i q 



The opposite orientation of the circular polarization is 
employed by substituting Cl — > — Cl in this definition. The 
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FIG. 2: The quasienergy spectrum for U = (top line) and 
U — 30meV (lower two lines) in both valleys, for weak radia- 
tion (left) and strong radiation (right). The color of the line 
indicates the weight of the static (n = 0) component. Thin 
red lines show the unirradiated spectrum. 



natural parameter by which to measure the strength of 
coupling of the electrons to the field is x — -Sm-i where 
F is the field intensity, Vt is the frequency of the radia- 
tion, and v F is the Fermi velocity. If x > 1 we say we 
are in the strongly irradiated regime. We take the dipole 
approximation and assume that the graphene is clean 
enough that we can ignore inter-valley scattering caused 
by defects such as lattice imperfections. We also neglect 
electron-electron interactions. The time dependent part 
of the Hamiltonian is periodic with period Tq — 2n/£l, 
so we can employ Floquet's theorem [14] to write the 
electron wave functions 'F(i) in the irradiated region as 



t ^ A {t) where the coefficient e A is the en- 



ergy of the dynamical state (called the 'quasienergy') 
The wave function in the temporal Brillouin zone <f> A (t) , 
defined for — 7r/f2 < t < tt/Q, is periodic in time and can 
be expanded over its Fourier components n and the state 
basis consisting of eigenfunctions of the static Hamilto- 
nian denoted tf> a . We therefore write 



l(x,*)= £ 



(i) 



We solve the time dependent Schrodinger equation for TL 
by taking the Fourier transform of the matrix elements 
of the Hamiltonian over the states i\) a and constructing 
the Floquet matrix [14]. Diagonalizing this matrix yields 
the quasienergies and the wave function coefficients Xna ■ 
In Figure 2 we show the low energy spectrum of irra- 
diated bilayer graphene with and without a static gap in 
each of the two valleys. The color of the line indicates 
the weight of the static component of the wave function, 
which represents the physically observable part of the 
dynamical state. In the left-hand column, the coupling 
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FIG. 3: (i) The total current, and (ii) the valley polarization 
as a function of the field parameters. In both plots, U = 
20meV, fj, = 12meV, and 77 = 0.3 x HQ, with fi = 2THz. The 
white contours denote the region of high valley polarization 
and significant current flow. 



parameter is x — 0.96 (weakly irradiated) while in the 
right-hand column x — 4.82 (strongly irradiated). We su- 
perimpose the unirradiated (F — 0) spectrum (red lines) 
for comparison. The radiation opens dynamical gaps at 
hfl/2 intervals (as was shown in the monolayer case [11]). 
Secondly, when there is a gate potential applied, dynam- 
ical states are present in the gapped region (see the lower 
two rows), and the quadratic shape of the low momen- 
tum part of the bands is restored for strong radiation. 
However, because K electrons couple more strongly to 
the radiation than K' electrons (due to the different sub- 
lattice composition of the wave functions), the weights of 
the static component of the Floquet states are drastically 
different in each valley. In the strongly irradiated regime, 
the notion of the static gap loses its meaning as there are 
many dynamical states with significant static component 
in that energy range. It is the dynamical states in the 
static gap which we utilize in the proposal for the valley 
filtering device. Reversing the polarization of the light 
or the orientation of U causes the K' valley to couple 
strongly. 

We now demonstrate the generation of valley polarized 
current by using irradiated bilayer graphene as a filter for 
an arbitrary current. We employ a tunnelling approach 
[15] where we suppose that the system consists of three 
parts, as shown in Fig. 1. They are the two graphitic 
'leads' described by Hamiltonians Hl,Hr = Hq with 
energy spectrum E a and chemical potential ±/x/2, and 
the central, irradiated region described by the time de- 
pendent Hamiltonian Hq — TC discussed above, with 
quasienergy spectrum ea and chemical potential fixed at 
zero. The contacts shown in Fig. 1 connect the graphene 
flake with external systems, and we do not consider their 
influence. The central region is linked to the leads via 
the coupling Hamiltonians Hql, Hcr- Denoting the op- 
erators for electrons in the leads by Ck Q i for i £ {L,i?}, 
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and the central area by dqA, we have 

H Ci = ]T ^ka,qAcLV + H - C - 

We assume that the central region is wide enough to for- 
bid electrons from tunnelling directly between the two 
leads. Since it has been shown [16] that transmission 
from bilayer graphene into gapped bilayer graphene is 
high for a wide range of the electron's angle of mo- 
mentum, we assume that for the transfer to occur, the 
electron's momentum is conserved and the energy of 
the states in the two regions must be sufficiently close. 
We parameterize this closeness by writing the function 
A(E) such that A(E) = for \E\ > r) and A(0) = 1 
so that T) describes the width of the allowed transition. 
Then, the coupling parameter is Vk a ,qA — V(5k, q A(£'ka — 
z^a) \Xoa\ ■ ^he quantity V has units of energy and pa- 
rameterizes the maximal strength of the coupling and we 
preserve the electron momentum via the S function. 

The valley component of the charge current in the 
right-hand lead is = — (dN^/dt), where is the 
number operator for the appropriate electron species. Us- 
ing a nonequilibrium Green's function analysis and tak- 
ing the steady state limit, we find that the current is 

^ = -|/||i>{ r «< 3<5 'M^ < 2 > 

where 3G r is the imaginary part of the full retarded 
Green's function in the central region, T contains the cou- 
pling parameters, and T = f c (E a ) — fn{E a ) depends 



on the distribution functions in the right lead and central 
region. The central region Green's function is calculated 
using the Floquet states derived above, and includes the 
self energy due to the two leads. 

To characterize the degree of valley polarization of the 
current, we define V = (J K — J K >)/{J K + J K ,) so that 
V = — 1(+1) corresponds to fully K' (K) polarized cur- 
rent. In Fig. 3 we plot the total current and the po- 
larization as a function of the radiation intensity and 
frequency for U = 20mcV and the chemical potentials of 
the leads arranged to drive current in the energy range 
corresponding to the static gap (// = 12meV). The area 
enclosed by the white contour shows where J > 0.04pA 
and V > 0.98 simultaneously, i.e. the region where the 
system parameters are tuned for significant current and 
very high polarization. Reversing the sign of U or the ori- 
entation of the polarization of the radiation leaves Fig. 
3(i) unchanged, but inverts Fig. 3(ii) so that the region 
of high current and polarization is in the K' valley. Iden- 
tification of the valley into which the current is polarized 
may be achieved by application of an in-plane electric 
field [7] which produces a valley-dependent Hall current 
which will result in a measurable asymmetry in the elec- 
tron density across the conducting channel. 

In summary, we have described the measurable char- 
acteristics of a graphene-based valley polarized current 
generator, where we expect current of ~ O.lpA and valley 
polarization of > 99%. Our work should provide neces- 
sary stimulus in the quest for valley tronics with graphene. 
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